
c --------------------------------------------------------------
c      program kpd1
c --------------------------------------------------------------
      parameter(Max_VMPTS=20000)
      parameter(Max_NODES=400)
      parameter(MAX_NBODY=2)
      
      complex z,v,vm, hub,z0, wall, zcr
      real*4 CHARD
      common /vortex/ NVORT,Z(Max_VMPTS),V(Max_VMPTS),VM(Max_VMPTS),
     >      gamma(Max_VMPTS)
      common /solid/ NBDIES,NWALL(2),WALL(Max_NODES,1),
     >   ZCR(Max_NODES,1), THETA(Max_NODES,1),
     >   HUB(2),z0(2),ninc(2),INC(15,2),XMAX(2)
      common/solid2/ CHARD
      common/system/ ndim,
     >   a(Max_NODES,Max_NODES),x(Max_NODES),ipvt(Max_NODES)
c --------------------------------------------------------------
c      gmtry
c --------------------------------------------------------------
ccc      common/system/ndim, a(Max_NODES,Max_NODES),x(Max_NODES),ipvt(Max_NODES)
c --------------------------------------------------------------
c       
c --------------------------------------------------------------
c  common /vortex/
C NVORT    NUMBER OF VORTEX PARTICLES
c z(*),v(*),vm(*),gamma(*)   data of vortex particle
c --------------------------------------------------------------
c  common /solid/
C NBDIES   NUMBER OF BODIES.
C NWALL    NUMBER OF WALL POINTS ON EACH OF THEM.
C WALL     ARRAY OF WALL POINTS. 
C ZCR      ARRAY OF CREATION POINTS.
C THETA    POLAR ANGLE OF ZCR POINTS.
C          THETA IS USED TO FIND IF VORTICES ARE
C          INSIDE SOLID. 
C ZO       USED ALONG WITH THETA.
C HUB(L)   IS THE HUB OF THE BODY "L". 
C INC      WILL HELP FIND OVER WHICH WALL POINT THE VORTEX IS,
C BY BISECTION. THE FIRST DIMENSION
C OF INC MUST BE AT LEAST LOG2(NWALL).
c --------------------------------------------------------------
c  common/system/
c    ndim
c    a(Max_NODES,Max_NODES)
c    x(Max_NODES)
c    ipvt(Max_NODES)
c
C A        IS THE MATRIX OF INFLUENCE COEFFICIENTS FROM
C          CREATION POINTS TO WALL POINTS.
c
c --------------------------------------------------------------
      integer nsysp
c      common/constpar/ nsysp
      parameter (nsysp =125)
c
c  common/constpar/
c       nsysp ---  系统矩阵最大容量
      
